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EXPERIMENTS  ON  HEAT  TRANSFER  OF  A  SUPERSONIC  GAS  FLOW  IN 


A  CIRCULAR  TUBE  AT  HIGH  TEMPERATURE  GRADIENTS 
V.  A.  Mukhin,  A.  S.  Sukomel,  and  V.  I.  Velichko 

An  experimental  investigation  was  made  of  heat 
transfer  during  the  flow  of  a  compressible  gas  in 
a  cooled  circular  tube  at  Mach  number  values  from 
0.1  to  4  and  T„  ./T„  from  1  to  3.1.  A  generalize- 

tion  of  the  experimental  data  Is  given. 

Most  studies  devoted  to  heat  transfer  in  a  flow  of  a  compressible 
gas  in  tubes  pertain  to  the  case  of  gas  motion  at  subsonic  velocities. 
Only  in  a  few  investigations  were  heat  transfer  and  drag  studied  dur¬ 
ing  the  flow  of  a  supersonic  gas  [4—6]  and  in  most  of  these,  the  meas¬ 
urements  were  carried  out  at  small  temperature  gradients.  We  made  an 
experimental  Investigation  of  heat  transfer  during  supersonic  air  flow 
in  a  circular  tube  at  large  temperature  gradients. 

Local  heat  transfer  was  measured  by  the  thick-walled  tube  method. 
The  experimental  heat  exchanger  was  a  thick-walled  brass  tube  575  mm 
long,  20.2  mm  in  inside  diameter  and  87  mm  in  outside  diameter.  The 
air  flowed  through  the  tube,  the  outside  of  the  tube  was  bathed  with 
water,  which  moved  through  the  annular  gap  between  the  casing  of  the 
heat  exchanger  and  the  outer  surface  cf  the  tube. 
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The  temperature  close  to  the  inner  surface  of  the  tube  was  meas¬ 
ured  at  20  points,  and  the  temperature  of  the  outer  surface,  at  6 
points.  Thermocouples,  which  were  used  to  measure  the  temperature 
close  to  the  inner  surface,  were  placed  in  radial  drillings  1. 5  mm  in 
diameter.  The  beads  of  the  thermocouples  were • soldered  to  the  bottom 
of  the  drilling  by  condensor  welding.  Openings  with  a  diameter  of 
0.25  mm  were  made  at  17  points  along  the  length  in  order  to  sample 
the  static  pressure. 

The  air  entered  the  experimental  section  from  the  stagnation 
chamber  through  interchangeable  nozzles.  The  experiments  were  carried 
out  with  five  nozzles,  of  which  one  was  subsonic  and  four  were  super¬ 
sonic,  calculated  for  Mach  number  values  of  2.  5#  5* 0,  5.5,  and  4.0. 

After  preliminary  dust,  oil,  and  moisture  removal,  the  air  passed 
through  a  noimal  orifice  which  served  to  determine  the  discharge. 
Further,  the  air  entered  an  electric  shaft  furnace  where  it  was  heated 
to  the  assigned  temperature  (200-800°C),  and  then  flowed  through  the 
chamber  in  which  the  flow  stagnation  parameters  were  measured.  A 
heater  was  wound  on  the  wall  of  the  stagnation  chamber  and  a  grid 
mixer  was  installed  at  the  chamber  entrance.  This  enabled  us  to  deter¬ 
mine  the  stagnation  temperature  with  a  sufficient  degree  of  accuracy. 
The  stagnation  temperature  was  measured  by  a  shielded  chromel-alumel 
thermocouple  and  the  stagnation  pressure  by  a  class  0.55  standard 
manometer. 

The  local  density  of  the  heat  flux  at  the  tube  wall  (excluding  a 
small  area  near  the  entrance  about  5  diameters  long)  was  determined 
under  the  assumption  that  the  temperature  field  in  the  wall  was  one- 
dimensional,  i.e.,  by  the  formula 

<?c  ='/lA  /. , 

_ i _  (1) 

Ii;  A’j/Ab 
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The  heat  conduction  of  the  wall  was  determined  experimentally 
for  each  thermocouple. 

The  flow  parameters  in  the  given  section  (velocity  coefficient  X 
and  flow  temperature  T)  were  found  on  the  basis  of  a  one-dimensional 
flow  model.'  The  stagnation  temperature  To  was  determined  from  the 
known  stagnation  temperature  at  the  entrance  and  from  the  distribution 
of  the  heat  fluxes  in  the  wall  from  the  entrance  to  the  section  under 
consideration.  The  heat-transfer  coefficient  was  calculated  by  the 
formula 


a  = 


(2) 


where 


r  =.  T 

J  n«c  —  4  0 


fc—  i 
'Hi 


(1  —  r)  X- 


/■  =  V  Pr®  —  7,16  - 10": 


Re°A 


7,“Fj"  1  *>r  *§—04-15; 


(3) 


f 


(  x 


1  4-0,0413 


X 


15  for  -p-  —  15427 ; 


\D)  '  '  \  D 

T*  =  T  +  0,5  (Te  -T)+  0,22  (7\,.c  -  T)  . 


Formula  (3 )  was  derived  from  the  generalization  of  the  experi¬ 
mental  data  [3]-  All  physical  properties  of  the  air  were  determined 
from  the  thermodynamic  temperature  of  the  flow  at  the  given  section. 

The  experiments  were  carried  out  at  Reynolds  number  values  from 
0.4  •  105  to  7.  O  '  105  and  Mach  number  values  at  the  tube  entrance 
from  0.1  to  4.0,  the  entrance  temperature  varied  from  200  to  800°C. 

An  analysis  of  the  experimental  data  showed  that  the  velocity 
and  pressure  distribution  over  the  tube  length  mainly  corresponded  to 
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the  one -dimensional  flow  model.  In  the  experiments  with  supersonic 
velocities,  we  observed  oblique  shocks  which  led  to  a  certain  disturb¬ 
ance  in  the  smooth  change  of  velocity  and  pressure  over  the  tube 
length.  The  magnitude  of  the  oblique  shocks  mainly  depended  on  how 
well  the  nozzle  was  coupled  with  the  experimental  tube.  Therefore, 
we  took  measures  to  couple  these  units  as  best  as  possible. 

The  character  of  the  lengthwise  change  of  the  local  heat-transfer 
coefficients  showed  that  under  certain  conditions  there  is  an  area  at 
the  start  of  the  tube  which  is  occupied  by  laminar  and  transition 
boundary  layers.  We  can  trace  that  the  length  of  this  area  is  re¬ 
tained  as  the  Re  number  increases.  Hereafter  we  will  examine  only 
points  corresponding  to  a  developed  turbulent  motion  and  to  complete 
thermal  stabilization  of  the  flow. 

At  first  we  carried  out  a  series  of  experiments  with  low  subsonic 
velocities  at  low  air  temperatures  at  the  inlet  (200-300°C). 

The  results  of  this  series  of  experiments  agree,  within  an  ac¬ 
curacy  to  5$  with  the  formula  of  M.  A.  Mikheyev  for  heat  transfer  in 
turbulent  motion  of  an  incompressible  fluid  in  tubes.  After  this,  we 
carried  out  experiments  at  high  temperatures  in  a  subsonic  nozzle  and 
experiments  on  supersonic  nozzles  both  at  low  and  high  temperatures. 

Figure  1,  which  was  calculated  by  means  of  parameter  t,  Ulus-  ' 
trates  the  effecv  of  gas  compressibility  on  heat  transfer: 


It  was  assumed  here  that  the  dependence  of  heat  transfer  on  the 
Reynolds  and  Prandtl  numbers  is  the  same  as  for  fluid  Incompressibility: 
Nu  ~  Re°»aPr0'43. 

As  we  see  from  Fig.  1,  the  Nusselt  number  varies  in  proportion 
to  t°>42.  Thus  the  dependence  of  the  Nusselt  number  on  t  was  somewhat 


greater  than  in  Petukhov's  study  [6],  which  led  to  a  difference  in 
Nusselt  number  values  which  did  not  exceed  10$.  This  small  divergence 
can  be  explained  by  the  different  experimental  conditions  in  these 
studies  (different  pressure  gradients,  temperature  gradients,  exponent 
In  the  temperature  dependence  of  viscosity). 

Figure  2  confirms  the  usual  dependence  of  heat  transfer  on  the 
Reynolds  number. 

As  a  result  of  treating  the  experimental  data,  we  obtained  the 
dependence 

Nu  =  O,022Reo-sPro'4:l  .  ( 4 ) 


The  formula  is  valid  for  cooling  of  a  turbulent  gas  flow  in  a 
tube  in  the  section  having  thermal  stabilization  at  values  of 

Re  =  (0,4 -r-7,0)  -  10s,  M  =  0,1  -h  3, 0*\  0  =  £us.  =  1  -f-3,1 . 

'  C 

The  scatter  of  the  experimental  points  relative  to  the  curve  cor¬ 
responding  to  formula  (4)  does  not  exceed  +10$,  which  is  within  the 
accuracy  of  the  experimental  data. 


Thus,  the  results  of  the  investigation  show  that  within  the  ac¬ 
curacy  of  the  experiments,  the  Nusselt  number  does  not  depend  on  the 
temperature  factor  in  spite  of  the  fact  that  the  latter  varies  in  wide 
limits  (from  1  to  j5. 1).  Of  course,  this  conclusion  is  valid  only  if 
the  gas  flow  is  cooled,  and  is  in  conformity  with  the  results  obtained 
in  earlier  studies  [1,  7]. 
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Fig.  2.  Dependence  of  the  Reynolds  (des¬ 
ignations  are  the  same  as  in  Fig.  l), 

B  =  Nu  . 

pr0,43TO,42 

DESIGNATIONS 

Ate  is  the  local  temperature  drop  at  the  tube  wall;  A  is  the  wall 
heat  conduction;  R0,  Ri,  R2  are  the  radius  of  the  tube  and  radii  of 
the  hot  Junction  inserts  of  the  inside  and  outside  thermocouples;  Tc 
is  the  wall  temperature;  T_  is  the  adiabatic  wall  temperature;  r 

at  v* 

is  the  temperature  recovery  coefficient;  Pr*  is  the  Prandtl  number 
value  taken  at  a  certain  temperature  T*  calculated  by  Eckert's  method. 
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SOLUTION  OF  AN  EQUATION  OF  GAS  MOTION  IN  A  RECTANGULAR 
CHANNEL  WITH  FRICTION  AND  HEAT  TRANSFER 

Yu.  I.  Danilov  and  B.  M.  Galitseyskiy 


A  method  is  given  for  the  solution  of  a  differential 
equation  of  gas  motion  in  a  rectangular  channel  under 
the  combined  influence  of  friction  and  heat  transfer. 


As  a  gas  moves  along  a  channel,  the  gas  parameters  at  a  given 
section  depend  on  the  friction  factor,  on  the  starting  parameters  of 
the  gas  in  front  of  the  channel,  and  also  on  the  law  and  intensity  of 
the  heat  effect. 

We  will  assume  that  the  law  and  the  form  of  the  heat  effect  are 
already  given,  therefore  the  stagnation  temperature  distribution  over 
the  length  of  the  channel  is  easily  determined  from  the  equation  of 
thermal  equilibrium 


dq  =  '/wFd(CpT). 


(1) 


The  equation  of  a  uniform  steady-state  gas  flow  in  a  rectangular 
channel  in  the  presence  of  friction  and  heat  transfer  has  the  follow¬ 
ing  form: 


d_ 

dx 


[P-f  psrJ+5 


o  wz 

2 zr 


=  0. 


(2) 
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Let  us  change  to  the  parameters  of  stagnation  and  the  gasdynamic 


functions 


P  -j-  p  iCI2  =  p  caz ). 


F' 


Then  Eq.  (2)  is  written: 

~  [*(>-)  VT}+±-±-t±-iVT  =  Q. 
dx  2  k  +  1  D 

To  solve  Eq.  (4)  we  will  make  a  change  of  variables 

s-i-  J^x/r  =  -A, 

DHl  ^ 

where  u  is  the  unknown  function  subject  to  determination. 
Then  Eq.  (4)  is  transformed: 


\ 


m(u  -j-  C)  ~  -j-  =  0, 

dx 


m 


J _ 

£>  !•  +  1  ' 


(3) 

(*) 

(5) 


(6) 


The  constant  of  integration  £  can  be  assumed  equal  to  zero, 
since  the  change  of  variable  r  =  v  -  c  reduces  Eq.  (6)  to  the  form 


dv_y- 
dx 


mv  ■ 


dv 

dx 


mT  =  0. 


Equation  (6)  at  a  constant  value  of  m  and  a  linear  change  of  the 
stagnation  temperature  over  the  channel  length  is  easily  integrated. 
Actually,  let 


r„  __  t, 

r  =  r  -i-  — — -i-  v 

1  'I'  ,  ’ 


(7) 


where  Ti  and  T2  are  respectively  the  stagnation  temperature  at  the 
channel  inlet  and  exit.  Then  Eq.  (6)  is  written: 


This  equation  will  be  integrated  by  introducing  a  parameter.  Let  us 
designate  =  z,  then 


,  /  t  _  r 

z~  -j-  muz  tir  ^7\  -} — z—j — — 


x  =  0. 


(9) 


Differentiating  Eq.  (9)  in  terms  of  u,  we  derive  after  appropriate 
transformations  a  first-order  linear  equation 

zu  ,  2z- 


Ull 

Hz 


,  ,  A  -  A 

z-  4-  m  — - - 

/ 


in 


(i!+ 


III 


A  -  A 
L 


=  0. 


(10) 


The  solution  of  this  equation  is  written  in  the  form 


u  = 


m 


~r 


A-A  ^-l  /  T„  —  A  \"  *' 

— - - Air-h 1  z  j  //t  — - - 


\ 


- A  \- ‘.vi  /  r.-Ai-*/. 

r)  - 


/ ,  r„  —  7\  \-  ‘/a 

■  c  1 2'  -r  /?:  — - L 1 

\  / 


(11) 


Eliminating  u  from  equations  (9)  and  (11),  we  obtain  an  equation 
relative  to  the  parameter  z: 


/«r 

z 


r  a 


z-  4-  m 


h=J\yh 


—  c 


z- 


*r  ni 


(12) 


We  will  determine  the  constant  of  integration  £  from  the  condition 
z  =  zi  when  T  •=■  Tx,  then  Eq.  (12)  is  transformed  to 


mT 


mTx 


\/ * 

/ 


-f  m 


A- A  •  A- A 


sir* 


I  T 
zj^m  - 


^r- 


z?  4- 


A- A 


+  r- -A.,^-1  z 


(1?) 


Since  the  parameter  £  is  uniquely  associated  with  X  and  T: 
z  »  *  -mx  then  substituting  its  value  into  Eq.  (ij5),  we  obtain 

the  following  expression: 


4-siah"1  /.,  A, 


-f-ciEh"*  X.  A., , 


(14) 
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where 


k+l  D\  Tx 


k~-  1  D\  T2 


or,  designating 


A*- 1- 


X2 


4-alnJr1  XA  =  0(X,  A), 


finally  for  determining  the  reduced  velocity  in  the  channel  we  will 
have  (Pig.  l) 


<!>(>•!.  A)  =  ,  A.). 


(16) 


Analogously  for  the  case  of  gas  cooling  according  to  the  linear  law, 
the  function 


CVA  1) 


i/¥ 


A1  -f  Arc  sine  A  A 


(17) 


Pig.  1.  Dependence  of 
$  on  X:  1,  2,  3,  4,  5, 
6)  Values  of  A  equal 
respectively  to  1,  0.9, 
0.8,  0.7,  0.6,  and  0.5* 


-11- 


Knowing  the  magnitude  of  the  reduced  velocity,  we  easily  deter¬ 
mine  all  remaining  gas  parameters  (static  temperature,  pressure,  den¬ 
sity)  by  using  the  gasdynamlc  functions.  For  example,  the  pressure 
losses  In  the  channel  will  equal: 

I  W- 1)  £/(*s)  J 

where 


r - 


When  the  change  of  stagnation  temperature  over  the  channel 
length  differs  from  the  linear  law  and  the  friction  factor  varies 
appreciably  of  the  channel  length,  we  recommend  the  method  of  divid¬ 
ing  the  channel  into  a  number  of  segments  In  which  we  can  assume  ap¬ 
proximately  a  constant  value  of  the  friction  factor  £  and  a  linear 
rise  of  the  stagnation  temperature.  Thus,  we  can  calculate  with  any 
degree  of  accuracy  a  channel  with  any  predetermined  law  of  heat  trans¬ 
fer. 

A 

For  sufficiently  small  values  X  «  1,  the  function  z(X)  «  75^  and 
Eq.  (6)  will  have  the  form 

mu  ~- +  m'T  =  0.  (18 ) 

Solving  thiB  equation,  we  derive  the  following  expressions  for 
determining  the  reduced  velocity  at  the  channel  exit  relative  to  the 
reduced  velocity  at  the  channel  inlet: 


(19) 


When  the  magnitude  X  »  1,  function  z(x)  can  be  assumed  approximately 

A  * 

equal  to  z(X)  Bs  ^  X.  Then  Eq.  (6)  is  transformed  as 

~-f  «ih-=0.  (20) 

ax 

It a  solution  is  written 

\f  jr  exp  -  J  nulx  .  (21 ) 

Formulas  (19)  and  (21)  are  derived  in  the  case  of  gas  flow  in  a 
channel  with  a  heat  conductor  and  respectively  reflect  the  character 
of  subsonic  and  supersonic  flow  wherein  when  we  derived  them  we  did 
not  permit  the  assumption  of  a  linear  change  of  stagnation  temperature 
over  the  channel  length,  i.e.,  they  are  valid  for  any  law  of  a  heat 
conductor.  However,  upon  using  the  obtained  relationships  in  practice, 
it  is  necessary  to  take  into  account  the  specific  conditions  of  the 
problem  and  the  required  accuracy  of  the  calculations. 

The  function  $  (X,  A)  depends  on  the  two  variables  A  and  X  and 
has  a  minimum  at  X  =  1.  Having  tabulated  the  function  $  (X,  A)  and 
compiled  a  table  for  a  wide  range  of  values  of  X  and  A,  we  can  fairly 
easily  and  rapidly  determine  the  distribution  of  the  reduced  velocity 
X  over  the  channel  length  and,  consequently,  all  other  parameters  for 
any  predetermined  law  of  heat  exchange. 

DESIGNATIONS 

^  is  the  heat  supplied  to  the  gas  in  the  channel;  F  is  the  flow 
area  of  the  channel;  C  is  the  heat  capacity  of  the  gas  at  constant 
pressure;  T  is  the  gas  stagnation  temperature;  p  is  the  density;  w  is 
velocity;  x  is  the  abscissa  readable  along  the  channel  axis  from  the 
channel  inlet;  P  is  static  pressure;  4  Is  the  friction  factor;  D  is 
the  hydraulic  diameter;  k  is  the  Index  of  the  adiabatic  curve;  X  is 
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the  reduced  velocity;  R  is  the  gas  constant;  c_  is  the  constant  of 
integration;  1  is  the  channel  length. 


HEAT  AND  MASS  TRANSFER  IN  EVAPORATION  PROCESSES 

A.  V.  Lykov 


The  effect  of" a  crossflow  of  a  substance  on  heat 
transfer  In  a  laminar  flow  past  a  moist  capillary- 
porous  plate  is  analyzed.  Regularities  in  heat  and 
mass  transfer  during  liquid  evaporation  from  capil¬ 
lary-porous  bodies  are  established.  It  Is  shown 
that  in  the  case  of  deepening  the  evaporation  sur¬ 
face,  the  heat-transfer  coefficients  are  larger 
than  those  in  evaporation  at  the  surface  of  a  body. 


Heat  and  mass  transfer  In  liquid  evaporation  from  an  exposed 
surface  and  from  capillary-porous  bodies  is  not  only  of  theoretical 
interest,  but  also  of  practical  value  for  engineering. 

Heat  and  mass  transfer  between  the  surface  of  a  liquid  and  the 
ambient  medium  (moist  air)  is  a  single  mutually  associated  process 
having  its  own  characteristics  which  differ  from  a  heat  transfer  proc¬ 
ess  complicated  by  mass  transfer.  An  even  more  complex  process  is 
heat  and  mass  transfer  of  a  moist  capillary-porous  body  with  the  flow 
of  a  heated  gas.  In  this  case,  heat  and  mass  transfer  of  the  body 
surface  with  the  ambient  medium  is  continuously  associated  with  heat 
and  mass  transfer  inside  the  capillary-porous  body  (the  interdepend¬ 
ence  of  the  external  and  internal  problems).  Such  processes  of  heat 
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and  mass  transfer  Include  those  of  drying  and  transpiration  cooling  by 
means  of  liquid  evaporation.  Heretofore  there  has  been  no  single 
opinion  not  only  on  the  physical  mechanism  of  transfer,  but  also  on 
the  qualitative  effect  of  mass  transfer  on  heat  transfer.  A  number 
of  investigators,  based  on  the  theory  of  gas  injection  into  the  bound¬ 
ary  layer  through  a  porous  plate,  consider  that  in  liquid  evaporation 
the  heat-transfer  coefficient  decreases  with  an  increase  in  evapora¬ 
tion  intensity,  whereas  in  condensation  of  vapor,  conversely,  the 
heat-transfer  coefficient  increases  with  an  increase  of  condensation 
intensity.*  In  the  first  case  the  crossflow  of  a  substance  is  di¬ 
rected  to  the  side  opposite  the  heat  flow  and  in  the  second  case  these 
flows  have  the  same  direction.  Then  the  boundary-layer  thickness  in 
evaporation  increases (the  boundary  layer  "swells"),  which  leads  to  a 
reduction  of  the  heat -transfer  coefficient.  During  condensation  the 
reverse  picture  occurs,  which  leads  to  an  increase  of  the  heat-trans- 
'fer  coefficient  with  an  increase  of  condensation  intensity. 

However,  the  experimental  investigations  of  A.  V.  Nesterenko  [l], 
G.  T.  Sergeyev  [2],  and  of  other  investigators  on  evaporation  of  var¬ 
ious  liquids  from  an  exposed  surface  showed  that  the  heat-transfer  co¬ 
efficients  in  evaporation  are  larger  than  those  without  mass  exchange 
(dry  heat  transfer)  under  the  same  hydrodynamic  conditions  and  temper¬ 
ature  differences.  This  difference  increases  with  an  increase  of  the 
relative  humidity  of  the  air.  When  drying  moist  materials,  the  heat- 
transfer  coefficients  are  larger  than  the  corresponding  coefficients 
of  a  dry  body  [3-6].  During  transpiration  cooling,  when  the  fluid  is 
continuously  fed  into  a  capillary-porous  body,  we  have  different  re¬ 
sults  for  different  investigators.  In  some  investigations  [7]  the  1 

*  Here  it  is  assumed  that  the  heat  needed  for  evaporation  Is 
mainly  transmitted  by  convection  from  the  heater  air. 
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heat -transfer  coefficients  increase  with  an  increase  of  the  evaporation 
intensity  and  in  others,  conversely,  they  decrease  [8,  The  results, 

obtained  experimentally,  are  of  considerable  interest.  In  our  article 
we  will  explain  them  on  the  basis  of  the  theory  of  liquid  evaporation 

from  capillary -porous  bodies. 

Transpiration  Cooling  in  Gas  Injection 
into  the  Boundary  Layer 


A  calculation  of  transpiration  cooling  by  the  method  of  boundary- 
layer  injection  through  porous  walls  was  studied  in  greatest  detail 
by  E.  Eckert  [10].  It  Is  based  on  a  solution  of  a  system  of  differ¬ 
ential  equations  of  heat  and  mass  transfer  for  a  laminar  boundary 
layer  with  flow  of  a  binary  gas  mixture  past  a  flat,  porous 
plate.  The  system  of  differential  equations  has  the  form: 


dn'x 

dx 


(1) 


d;ox  ,  dwx 
l >wx  -f  pwy  -z~- 
dx  dy 


dt  . 

Cp  P  “-'.V  '\j  -f  CP  P  WU 


dt 


-r  P  D(c„:  —  cpl) 


d  f. 

dwx  \ 

dy  V1 

dy  )’ 

5  ( 

.  dt 

dy  { 

dy 

dt 

dy 
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(2) 

(3) 


dx 


dp„, 

dy 


dy 


dy  ) 


(*) 


Thermal  diffusion  (Soret  effect)  and  the  diffusion  thermo  effect 
(Dufour  effect)  are  disregarded  as  small  magnitudes. 

The  following  boundary  conditions  were  taken: 

’•’ton  //  =  0  —  0,  wv  =  Wy,  t  =  tr  [>1#  =  p10»;,  (5 ) 
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(6) 


wl*n  !i  ”  '-0  wx  —  ^c<  l  —  ‘c<  Pm  —  Pioc- 


In  addition.  It  was  assumed  that  the  linear  transverse  velocity 
of  the  change  along  the  surface  wg  (In  the  x-direction)  Is  inversely 
proportional  to  «/x. 

The  results  of  the  calculations  for  the  case  where  the  mass  flow 
of  the  injected  gas  is  directed  from  the  body  surface  (analog  of  the 
evaporation  process)  are  shown  in  Fig.  1.  We  see  from  this  figure 
that  the  heat-  and  mass -transfer  coefficients  decrease  with  an  in¬ 
crease  of  the  parameter 


Z  = 


/Re*, 


(7) 


where  Ji  =  pwg  is  the  intensity  of  mass  transfer  in  a  direction  normal 
to  the  wall  surface  (the  intensity  of  the  mass  crossflow  of  a  sub¬ 
stance). 

Consequently  the  heat-transfer  coefficient  a  decreases  with  an 
Increase  of  Intensity  of  the  crossflow  of  the  substance.  The  ratio 
a/a0  =0.9  when  Z  =0.05  (Fig.  l).  Consequently,  when  Z  =  0.05  the 
decrease  In  the  heat-transfer  coefficient  Is  10$.  A  twofold  decrease 
(a/a0  =  0.5)  corresponds  to  a  value  of  Z  «  0.5.  In  the  study  by 
Shulman  [9]  with  transverse  flow  past  a  hollow,  porous  ceramic  cylin¬ 
der,  inside  which  water  was  fed  to  certain  regions,  the  parameter  Z 
varied  from  0. 015  to  0.05.  The  evaporation  intensity  was  from  5.6  to 
18.6  kg  •  m-2  •  hr-1  in  the  Reynolds  number  range  from  2  •  104  to 
8  •  104  and  with  an  air  temperature  change  from  JO  to  150°C.  Conse¬ 
quently,  a  decrease  in  the  heat-transfer  coefficient  caused  by  cross- 
flow  was  less  than  10$,  i.e.,  it  lies  within  the  accuracy  of  the  ex¬ 
periment.  Moreover,  a  decrease  In  the  local  heat-transfer  coefficient 
was  noted  upon  an  increase  in  the  evaporation  Intensity  by  a  factor  of 
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Fig.  1.  Dependence 
of  relative  heat- 
transfer  coefficient 
(a/a0)  on  parameter 
Z  based  on  Eckert's 
data  [10]:  l)  Pr  = 
=  Sc  =  1;  2)  Pr  = 


1. 5-1*9/  i.e.,  by  50-90$.  If  we  take  into  account  that  the  process 
of  liquid  evaporation  is  analogous  to  the  process  of  blowing  air  into 
the  boundary  layer,  then  under  conditions  of  adiabatic  evaporation 
the  parameter  Z  characterizing  the  effect  of  crossflow  of  a  substance 
on  heat  and  mass  transfer  can  be  presented  in  a  different  form.  The 
intensity  of  evaporation  is 

A  =  — —  (A  -  4)  =  —  Xu,  (/,  - 4).  (8 ) 

r  x 


In  a  laminar  flow  of  moist  air  past  a  plate  Nu„  -  k  x/Re.  where  k 

AAA  .A. 

is  a  constant.  Then  the  parameter 


A*  c„  1  / 

Pr  r 


(9) 


Consequently,  Z  does  not  depend  on  the  velocity  of  air  motion, 
but  depends  on  the  psychrometric  difference. 

Let  us  make  an  approximate  calculation.  For  Z  =  0. 05  the  temper¬ 
ature  difference  At  will  be  At  =  260°C.  In  the  calculation  we  assumed 
'tjj  «■*  30°C,  r  <■  579  kcal  •  kg-1,  Ax  *  0.53.  Consequently,  only  at  a 
temperature  difference  of  At  =  260°  will  the  reduction  of  the  heat- 
transfer  coefficient  be  of  the  order  of  10$. 

Similar  results  were  obtained  earlier  [8]  with  the  flow  of  hot 

air  past  a  porous  plate.  The  evaporation  intensity  in  these  experl- 

—2  -1 

ments  did  not  exceed  22.7  kg  •  m  •  hr  .  The  Reynolds  number  varied 
from  104  to  2  •  106,  and  the  air  temperature  from  18  to  140°.  Evapor¬ 
ation  was  under  adiabatic  conditions.  The  temperature  difference  did 
-not  exceed  79°*  Consequently,  crossflow  of  a  substance  cannot  render 
a  perceptible  effect  on  the  reduction  of  the  heat-transfer  coefficient. 
However,  the  authors  established  that  there  is  a  dependence  for  the 
average  Nusselt  number: 

Nu  =  0,00-155  Re0'8  ^ — — -JU.  (10) 


Therefore,  the  reduction  of  the  heat-transfer  coefficient  was 
more  than  threefold. 

Although  this  analysis  of  the  study  [8]  is  of  a  somewhat  arbi¬ 
trary  nature  since  evaporation  of  water  from  a  porous  plate  occurred 
in  a  turbulent  flow,  the  decrease  in  the  heat-transfer  coefficient 
with  an  increase  of  the  temperature  difference  At  is  the  result  of 
other  causes.  •  We  must  point  out  that  in  a  turbulent  flow,  the  effect 
of  crossflow  of  a  substance  on  the  change  of  the  heat-transfer  coef¬ 
ficient  will  be  less  than  in  a  laminar  flow. 

Drying  of  moist  materials  in  convective  driers  is  done  approx¬ 
imately  under  adiabatic  conditions.  The  temperature  differences  in 
hot-air  drying  is  considerably  lower  than  250°  in  the  overwhelming 
majority  of  cases.  Therefore,  crossflow  of  a  substance  through  the 
boundary  layer  produced  by  moisture  evaporation  has  virtually  no  ef¬ 
fect  on  the  value  of  the  heat-transfer  coefficient.  Numerous  experi¬ 
ments  established  that  at  a  constant  drying  rate  the  heat-transfer 
coefficient  is  considerably  larger  than  in  heat  transfer  of  a  dry 
body,  all  other  conditions  being  equal.  Starting  at  the  critical 
moisture  content,  the  heat-transfer  coefficient  decreases  during  the 
course  of  drying  time,  gradually  approaching  the  value  of  the  heat- 
transfer  coefficient  of  a  dry  body.  Consequently,  the  heat-transfer 
coefficient  decreases  with  a  decrease  in  the  drying  rate. 

Liquid  Evaporation  from  an  Exposed  Surface 

In  liquid  evaporation  from  an  exposed  surface,  the  heat -transfer 
coefficient  increases  with  an  increase  of  evaporation  intensity.  This 
effect  of  the  evaporation  process  on  heat  and  mass  transfer  is  char¬ 
acterized  by  the  Gukhman  number  which  is  the  thermodynamic  criterion 
of  evaporation. 
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The  calculation  formulas  in  a  criterial  form  are: 


Ku  =  /IPr^RcfGu"1, 

(11) 

Mi'  =  /■i'Sc°-33Rc"'Gu"1', 

(12) 

where  A,  A',  n,  n1,  m,  m'  are  constants,  determinable  by  experiment, 
dependent  on  the  hydrodynamic  flow  conditions  (Reynolds  number  range). 
For  example,  according  to  the  data  of  G.  T.  Sergeyev  [2],  in  the 
Reynolds  number  range  from  4  •  104  to  16  •  104,  the  constants  in 

formulas  (11)  and  (12)  are  equal  to:  A  =  0.0 86;  A'  =  0.094;  n  =  n1  = 

=  0. 8;  m  =  m'  =  0.  2. 

The  most  probable  explanation  of  the  evaporation  effect  on  the 
intensity  of  heat  and  mass  transfer  is  the  hypothesis  of  volume  evap¬ 
oration  associated  with  the  dynamic  character  of  the  processes  of 
sorption  and  desorption. 

The  essence  of  this  hypothesis  is  that  minute  liquid  droplets 
impinge  upon  the  boundary  layer.  The  main  cause  of  the  separation  of 

droplets  from  an  exposed  surface  of  a  liquid  is  the  presence  of  the 

processes  of  focal  condensation  and  the  interaction  of  the  gas  flow 
with  the  liquid  surface.  According  to  de  Bour's  dynamic  theory  of  ad¬ 
sorption  [11],  the  evaporation  process  is  a  dynamic  process  of  desorp¬ 
tion  and  sorption.  The  molecules  of  the  liquid  not  only  leave  the 
surface  (evaporation)  but  are  continuously  returned  (condensation). 

The  evaporation  rate  is  proportional  to  the  difference  of  the  molecu¬ 
lar  flows  leaving  and  returning  to  the  liquid  surface.  The  investiga¬ 
tions  of  N.  N.  Fedyakin  [12]  showed  that  condensation  does  not  occur 
uniformly  along  the  surface,  but  on  certain  sections,  and  incomplete 
wetting  of  the  liquid  surface  by  the  adsorbed  layer  of  liquified  vapor 
takes  place.  In  the  condensation  regions  drops  are  formed  which. 


-21- 


being  less  strongly  bound  with  the  liquid,  are  carried  away  by  the  air 
flow  to  the  boundary  layer. 


The  evaporation  of  drops  in  the  boundary  layer  we  will  call 
volume  evaporation,  it  is  a  volume  source  of  vapor  and  a  negative 
source  of  heat  in  boundary-layer  equations.  In  the  presence  of  vol¬ 
ume  evaporation,  the  right-hand  part  of  Eq.  (3)  will  contain  a 
third  term  equal  to  rl,  where  I  is  the  volume  capacity  of  the  vapor 
source  (kg  •  m-3  ■  hr).  The  right-hand  part  of  Eq.  (4)  must  be  sup¬ 
plemented  by  the  magnitude  of  the  positive  vapor  source  (I). 

By  methods  of  the  similarity  theory,  we  will  find  from  Eq.  (3) 
the  dimensionless  variable 


(15) 


The  magnitude  rl  equals  the  amount  of  heat  which  is  needed  for  volume 
evaporation.  If  we  designate  in  terms  of  Ny  the  number  of  drops  per 
unit  volume  of  the  boundary  layer,  and  is  the  mean  radius  of  a  drop, 
then  an  elementary  calculation  leads  to  the  relationship 


K,  =  .  (14) 

Tc- 

The  magnitude  8ttRNv12  depends  on  the  physical  properties  of  the  liquid 
and  the  hydrodynamics  of  the  flow. 

The  temperature  difference  At  =  [t(x,  y)  -  t  ]  varies  in  the  di¬ 
rection  of  the  x-  and  y-coordinates.  The  relative  value  At/Tc  char¬ 
acterizes  the  local  thermodynamic  intensity  of  evaporation.  The  max¬ 
imal  value  of  this  magnitude  equals  the  Gukhman  number  (Gu  = 

=  (Tq  -  Tm)/Tc).  From  this  point  of  view,  Gu  characterizes  the  poten¬ 
tial  possibility  of  humid  air  in  volume  evaporation.  The  hypothesis 
of  volume  evaporation  is  in  need  of  confirmation  by  direct  experiments. 
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The  second  cause  of  the  intensification  of  the  heat-transfer  process 
in  evaporation  is  the  disturbance  in  the  near-wall  boundary  layer  by 
focal  evaporation  processes.  In  liquid  evaporation  at  ordinary  temper¬ 
ature  and  pressure,  the  volume  of  the  substance  increases  by  about  103 
times,  and  on  condensation  of  the  vapor  the  same  decrease  in  volume 
occurs.  As  a  result  of  focal  processes  of  evaporation  and  condensation 
the  structure  of  the  laminar  boundary  layer  is  disturbed,  which  leads 
to  an  intensification  of  heat  and  mass  transfer.  This  effect  is  espe¬ 
cially  demonstrated  upon  evaporation  in  a  vacuum,  when  the  change  of 
volume  in  phase  transformations  reaches  an  order  of  106.  This  leads 
to  an  increase  in  the  heat -transfer  coefficient  by  about  one  order. 

The  hypothesis  of  stream  heat  and  mass  transfer  in  evaporation 
in  a  vacuum  was  developed  by  A.  A.  Gukhman  [15].  It  is  completely 
natural  that  these  effects  intensifying  heat  and  mass  transfer  take 
place  in  evaporation  under  conditions  of  ordinary  barometric  pressure 
although  their  influence  will  be  appreciably  smaller. 

When  drying  moist  materials  the  effect  of  volume  evaporation  on 
the  heat-  and  mass-transfer  process  will  be  smaller  as  compared  with 
liquid  evaporation  from  a  free  surface.  However,  the  Intensification 
of  heat  transfer  by  the  effect  of  focal  evaporation  will  evidently  be 
greater.  This  is  because  sorption  and  desorption  processes  take  place 
on  the  surface  of  macro-  and  microcapillaries  in  capillary-porous  bod¬ 
ies.  It  is  important  to  note  here  that  the  external  heat  and  mass 
transfer  of  moist  materials  Is  continuously  associated  with  physico¬ 
chemical  processes  on  the  surface  of  a  capillary-porous  body.  Deepen¬ 
ing  of  the  evaporation  surface  is  the  principal  factor  influencing  heat 
and  mass  transfer  between  moist  capillary-porous  bodies  and  hot  gas. 
Thus  heat  and  mass  transfer  of  a  body  surface  with  the  ambient  medium 
should  be  considered  as  a  combination  of  heat  and  mass  transfer  in 
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the  boundary  layer  of  humid  air  and  in  the  boundary  layer  of  a  capil¬ 
lary-porous  body  (evaporation  zone). 

Liquid  Evaporation  from  Capillary-Porous  Bodies 


To  analyze  the  complex  problem  of  heat  and  mas3  transfer  in  evapor¬ 
ative  transpiration  cooling,  we  will  use  0.  Krisher's  method  [14]. 

The  essence  of  this  method  is  as  follows. 

If  we  neglect  the  effect  of  a  crossflow  of  a  substance,  the  dif¬ 
ferential  equation  of  heat  transfer  in  the  boundary  layer  with  laminar 
flow  past  a  flat  plate  can  be  written  as*: 


dl  =  a  clH 
dx  dtf 


The  boundary  conditions  are 


(15) 


when  jl  —  0  /  (X,  0)  =  ls%  when  -V  =  0  /  (0,  //)  —  /c, 
wh-'-n  1J  — *  co  /  (.v,  co)  =  /c. 


(16) 


The  flow  velocity  wx  is  a  function  of  the  coordinates  which  is  deter¬ 
mined  from  the  solution  of  the  equation  of  motion.  Krisher's  method 
assumes  that  wx  is  a  constant  magnitude  and  equal  to  the  average  flow 
velocity  in  the  boundary  layer  (w  -  wx).  In  actual  processes  a  con¬ 
stant  velocity  occurs  only  with  liquid  flow  without  friction,  i.e., 
at  a  very  small  coefficient  of  internal  friction.  In  the  case  of  a 
viscous  liquid  such  an  assumption  (w  =  wx  const)  Is  a  method  for 
solving  the  problem  of  heat  transfer  In  the  boundary  layer. 

In  the  case  of  flow  about  an  infinitely  long  plate  (1  -*■  «°)  the 
solution  of  differential  equation  (15 )  with  boundary  conditions  (16) 
has  the  form 


& 


*  Heat  transfer  by  vapor  diffusion  can  be  disregarded  as  a  small 
magnitude. 
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=  erf 


(17) 


j  (-V.  U)  —  ‘  s 

t'-t* 


The  local  Nusselt  number  is 


f-^%- 

\  Wax 


X  (It  (a-,  0) 

(4-g 


(18) 


Differentiating  solution  (17)  with  respect  to  £  and  assuming  y  =  0, 
we  obtain 


The  average  Nusselt  number  over  the  surface  is 


/ 


(19) 


(20) 


Fig.  2.  Diagram  of 
the  calculation  of 
heat  transfer  in  * 
evaporative  trans¬ 
piration  cooling. 


In  order  to  compare  the  obtained  results  with  known  formulas  for 
the  Nusselt  number  in  a  laminar  flow  past  a  plate,  we  must  determine 
the  magnitude  wx-  If  the  velocity  profile  wx(y)  is  assumed  a  cubic 
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parabola,  the  mean  integral  velocity  [15]  is 

_  j  <> 

k'x  —  1  -',(7)  <(■•'  -’-a»c,  (21) 

j  8 

where  5  is  the  boundary-layer  thickness. 

Then  for  moist  air  (Pr  -  0.7) 

Nu  «  0,625  ™=  ]■ ■  'Pc  «  0.74  Re  .  (  22 ) 

This  result  differs  only  by  12$  from  the  known  empirical  formula 
Nu  =  0.66  */r 5".  The  assumption  of  the  constancy  of  velocity  wx  in 
the  solution  of  the  differential  equation  of  the  boundary  layer  is 

thus  fully  admissible. 

Our  problem  concerning  evaporative  transpiration  cooling  can  be 
set  up  as  in  Fig.  2. 

The  differential  equation  of  heat  transfer  stays  the  same  (15). 
The  boundary  conditions  will  be: 


t  (0,  y)  =  4,  t  (a-,  00).  =  4,  i(x,  •)  =  4.; 


(25) 


'  (0.  x)  _  ■  fa  (0.  -v) 

<)y  T  dy 


—  \l  (0,  .v)  -41. 


(24) 


Here  it  is  assumed  that  the  temperature  at  the  evaporation  surface 
equals  the  temperature  of  a  wet -bulb  thermometer.  Owing  to  the  small 
boundary-layer  thickness  of  a  body  the  temperature  distribution  in  it 
follows  the  linear  law.  In  this  case  (24)  we  can  write: 


_  +  tl  [i  (0,  x)  -  41  =  0.  ( 25 ) 

oil 

h 

where  H  =  r—  is  some  magnitude  analogous  to  the  relative  heat  transfer 
coefficient. 

The  solution  of  differential  equation  (25)  with  boundary  condi¬ 
tions  (23)  and  (25)  has  the  form 
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(26) 


/  (.V.  //)—  /M 


or  I 


A  ' 


2  !  «.v 


.  /  .'/  7  ~'.v 

eric  _ 

I  2  1 /  «.v 


Prom  (26)  we  derive  the  solution  of  (17)  as  a  particular  case.  If 
evaporation  occurs  at  the  body  surface  (£  =  0),  then  the  second  term 
of  the  right-hand  part  of  (26)  equals  zero,  because  when 

c-0 

The  temperature  at  the  body  surface  (y  =  0)  will  not  be  constant, 
but  varies  along  thfe  x-axis: 

exp  (7/- j  orfc  ^  j  /  -|LJ  .  (27) 


4  --  /» 


At  the  plate  edge  (x  =  0)  the  temperature  of  the  body  surface 
ts  [tg  =  t(x,  0)]  equals  the  air  temperature  (tg  =  tc),  whereas  at  an 
appreciable  distance  (x  -*•  °°)  it  equals  the  temperature  of  the  wet- 
bulb  thermometer  (tg  =  tm).  Consequently,  the  temperature  difference 
At  (At  =  t„  -  t„)  varies  from  zero  at  the  plate  edge  to  a  constant 

C  S  .  I  -  .  ,1  — 

magnitude  (tQ  -  tm).  This  is  a  very  important  fact  determining  the 
specific  characteristics  of  heat  and  mass  transfer  upon  deepening  of 
the  evaporation  surface  into  the  body.  If  evaporation  occurs  at  the 
body  surface,  then  under  adiabatic  conditions  the  body  surface  temper¬ 
ature  is  constant  and  equals  the  temperature  of  the  wet-bulb  thermom¬ 
eter. 

It  is  known  from  the  general  theory  of  heat  transfer  that  if  the 
temperature  difference  At  Increases  in  the  flow  direction,  the 
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heat-transfer  coefficient  becomes  larger  than  at  a  constant  temperature 
[14].  Consequently,  upon  deepening  the  evaporation  surface  the  heat- 
transfer  coefficient  is  larger  than  when  evaporation  occurs  at  the 
surface.  If  in  the  first  approximation  we  assume  that  the  heat- 
transfer  coefficient  in  evaporation  at  a  body  surface  equals  the  heat- 
transfer  coefficient  of  a  dry  body,  then  on  drying  with  a  deepening 
of  the  evaporation  surface  the  heat -transfer  coefficient  will  be  larger 
as  compared  with  the  heat-transfer  coefficient  of  a  dry  body.  This 
increase  of  the  heat-transfer  coefficient  must  be  reflected  in  the 
calculation  formulas  Nu  =  f(Re)  by  introducing  an  additional  general¬ 
ized  argument.  Since  At  is  the  cause  of  the  change  in  the  heat- 
transfer  coefficient,  then,  naturally,  the  parametric  criterion  (gen¬ 
eralized  argument)  will  be  the  Gukhman  number  or  •  We  will  dwell 

on  this  in  detail. 

The  local  Nusselt  number  is 


Nu,  = 


dl  (x,  0) 


/c  —  /(0,  x)  d,j 

X  (1  —  exp  A'-erfc  A'j~‘  . 


—  ]/" Pe,  A' exp  A'2 erfc /<  X 


(28) 


where  the  dimensionless  variable 


V  Pe: 


(29) 


characterizes  the  effect  of  deepening  the  evaporation  surface  on  heat 
and  mass  transfer  of  capillary-porous  bodies.' 

We  will  designate 


l\K)-=  j  -  A' exp  A"2  erfc  A'. 


(30) 


Then 


y  _  J..  L_'i> 

i 


~  Nu, 
Pe, 


/  (A) 


1  -  Wtm 


(31) 
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We  derive  (19)  from  the  solution  of  (31).  If  evaporation  occurs  at 
the  surface  (K  =  »),  then  N  =  1  since  f(K)  =  1,  i.e., 

i  “ 

which  fully  agrees  with  solution  (19). 

Thus  the  dimensionless  magnitude  N  characterizes  the  relative 
Increase  of  the  local  Nusselt  number  during  evaporation  of  moisture 
from  capillary-porous  bodies  as  compared  with  evaporation  of  moisture 
at  the  body  surface. 

Figure  3  shows  a  graph  of  N  =  f(K)  from  where  we  see  that  the 
magnitude  N  decreases  with  an  increase  of  K,  gradually  approaching 
unity. 

In  the  range  of  K  values  from  0.1  to  5  (0.1  <  K  <  5.0)  the  depend¬ 
ence  N  =  f(K)  can  be  approximated  by  the  relationship 

A'  —  1,31  A'”0,12.  (32) 

It  is  known  from  the  theory  of  drying  moist  materials  that  the 
distance  of  the  evaporation  surface  from  the  body  surface  in  the  con¬ 
stant  drying  rate  period  in  the  first  approximation  is  directly  pro¬ 
portional  to  the  psychrometric  difference  At (At  =  t  -  t  ).  Then  the 

c  in 

dimensionless  variable  K  will  be  inversely  proportional  to  (tc  -tm) 
and  consequently  inversely  proportional  to  the  Gukhman  number  (K  ~ 

~  Gu_1).  Hence  it  follows  that  N  ~  Gu0'1,  which  takes  place  In  exper¬ 
iments  on  heat  and  mass  transfer  in  the  drying  process  [2]. 

Only  at  the  constant ’drying  rate  period  is  the  temperature  of  the 
evaporation  surface  t(x,  -£)  constant,  starting  with  the  critical 
moisture  content  its  temperature  increases  with  time  of  drying,  grad¬ 
ually  approaching  the  air  temperature  which  it  achieves  at  equilibrium 
moisture  content.  Hence  it  follows  that  Nux  in  the  falling  drying 
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rate  period  will  decrease  with  time,  gradually  approaching  the  Nusselt 
number  for  a  dry  body. 


Fig.  3-  Dependence  of 
coefficient  N  on  param¬ 
eter  K. 


It  is  of  interest  to  determine  approximately  the  magnitude  £. 

For  a  porous  ceramic  (X  =  0.2  kcal  •  m_1  •  hr-1  •  deg-1)  when  Re  = 

=  6  •  104  and  for  x  =  20  mm  and  parameter  K  =1.5,  the  magnitude 
£  =  0. 5  mm.  In  this  case,  according  to  the  graph  in  Fig.  3,  coef¬ 
ficient  N  =  1.25,  i.e.,  the  heat-transfer  coefficient  in  drying  is 
approximately  25$  greater  than  that  for  a  dry  body,  if  we  consider 
that  heat  transfer  in  evaporation  at  the  surface  of  a  body  is  identi¬ 
cal  to  heat  transfer  of  a  dry  body.  For  the  parameter  K  =  0. 25,  the 
magnitude  £=  3  mm,  and  the  coefficient  N  =  1.55,  i.e.,  the  heat  trans¬ 
fer  coefficient  is  about  25$  greater  during  drying  as  compared  with 
heat  transfer  of  a  dry  body. 

It  is  completely  natural  that  at  small  values  of  £  it  is  prac¬ 
tically  impossible  to  measure  accurately  the  temperature  of  the  body 
surface.  Thermocouples  embedded  on  the  "surface"  of  a  body  for  all 
practical  purposes  show  the  wet-bulb  temperature.  Therefore  the  heat- 
transfer  coefficient  can  be  considered  as  the  ratio  of  heat  flow  to 
the  psychrometric  differences  (t  -  t_): 

C  Iu 
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Then  the  local  Nusselt  number  is 


A 


_ -V _ (kjx.  0) 

(I'-ij  ():/' 


(34) 


After  simple  transformations  we  obtain 


V  l  l  \- 


(35) 


where  Nm  is  the  coefficient  showing  the  relative  change  of  the  number 
Nu^  and  consequently,  of  the  coefficient  owing  to  deepening  of 
the  evaporation  surface.  The  graph  Nm  =  f(K),  given  in  Fig.  4,  shows 
that  the  coefficient  N  increases  with  increase  of  K. 


tV  cs  i"  t.  6  r 


Fig.  4.  Dependence  of  the  coef¬ 
ficient  Nm  on  parameter  K  in  or¬ 
dinary  and  logarithmic  scales. 


Since  the  dimensionless  variable  K  is  inversely  proportional  to 
the  psychrometric  difference  (tc  -  t  ),  the  Nusselt  number  will  de¬ 
crease  with  its  increase  or  that  'of  the  Gukhman  number. 

In  a  small  variation  range  of  K  the  dependence  Nm  =  f(K)  can  be 


presented  as 


N  =  BKb, 
m 


(56) 


where  B  and  b  are  constants  (0  <  q  <  l).  For  example,  in  the  range  0. 3  < 
<  K  <  1.  5  B  =  0.73,  b  -  0.46,  and  when  1.  5  <  K  <  5  B  -  0.80,  b  -0.15 
(Fig.  4). 

o 

Consequently,  in  the  range  0.  3  <  K  <  1.  5  the  coefficient  Nm  will 
be  directly  proportional  to  Gu”0,46.  A  similar  relation  occurred  in 
Shull 'man's  work  [9]  in  which  he  experimentally  obtained  the  relation¬ 
ship 


/Re, 


=  /'(  X)Gu-o-*, 


(57) 


where  f(X)  is  a  function  of  a  dimensionless  coordinate. 

By  its  structure  Nm  =  0.73K0,46  is  close  to  empirical  formula 
(10),  since  the  variable  K  is  inversely  proportional  to  At.  We  can 
note  that  the  authors  of  article  [8]  carried  out  their  experiments 
using  the  Gukhman  number.  In  this  treatment  they  obtained  the  fol¬ 
lowing  dependence: 


N'n  =  0.00695  Re"-'1  G11  ,  (38  ) 

which  almost  agrees  with  the  calculation  formula  with  respect  to  the 
exponent  £  in  spite  of  the  different  conditions  of  flow  past  the 
plate.  We  should  point  out  that  a  comparison  of  these  formulas  Is 
tentative  since  the  magnitude  £  depends  on  the  capillary-porous  struc¬ 
ture  of  the  body,  its  physicochemical  properties,  and  in  the  general 
case  is  a  function  of  the  parametric  criterion  (Tc/T  ).  However 
formulas  (3i),  (55)  convincingly  show  that  in  these  studies  [8-9]^ 
evaporation  of  water  occurred  at  a  certain  depth  from  the  body  surface 
and  the  decrease  in  the  heat-transfer  coefficient  with  an  increase  of 
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evaporation  rate,  established  from  experimental  data,  is  explained 
not  by  the  crossflow  of  mass  to  the  boundary  layer  but  by  the  method 
of  calculating  the  value  of  the  heat-transfer  coefficient.  The  heat- 
transfer  coefficient  in  drying  of  capillary -porous  bodies  or  in  evapor¬ 
ative  transpiration  cooling  is  greater  than  the  heat-transfer  coef¬ 
ficient  of  a  dry  body  at  ordinary  temperature  differences  (At  =  • 

=  250-300°C). 

Actually  in  these  studies  [8-9]  not  the  heat-transfer  coefficient 
ax  was  determined  but  the  over-all  heat  transfer  coefficient 

(39) 

\  »•  '•  / 

It  is  completely  natural  that  the  coefficient  k  decreases  with  an  in¬ 
crease  of  the  psychrometric  difference  (t  -  t  )  since  £  in  the  first 
approximation  is  proportional  to  (tQ  -  tm). 

The  basic  conclusion  of  our  investigation  is  that  heat  and  mass 
transfer  of  capillary-porous  bodies  with  the  ambient  medium  Is  a 
single  interrelated  process  of  heat  and  mass  transfer  in  the  boundary 
layer  of  the  body  and  In  the  boundary  layer  of  the  medium. 

CONCLUSIONS 

In  this  article  we  have  analyzed  Eckert's  solution  [10]  on  the 
effect  of  mass  transfer  on  heat  transfer  with  gas  Injection  Into  a 
laminar  boundary  layer.  It  was  shown  that  with  liquid  evaporation 
from  a  porous  plate,  a  decrease  In  the  heat  transfer  coefficient  of 
more  than  10$  can  take  place  only  at  a  temperature  head  At  j>  260°C. 

Under  the  usual  conditions  of  drying  and  evaporative  transpira¬ 
tion  cooling  this  effect  can  be  neglected.  An  analysis  of  the  exper¬ 
imental  data  on  liquid  evaporation  from  a  free  surface  shows  that  an 
additional  argument,  the  Gukhman  number,  should  be  Introduced  into  the 
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impirical  relations  Nu  =  f(Re,  Pr)  (see  formulas  (11),  (12)).  The 
Gukhman  number  characterizes  the  thennodynamic  intensity  of  volume 
evaporation  in  a  boundary  layer. 

The  problem  of  heat  transfer  in  evaporative  transpiration  cooling 
with  the  presence  of  deepening  of  the  evaporation  surface  is  solved 
by  Krlsher's  method  [14]  (see  Eq.  (26)).  It  was  shown  that  in  this 
case  the  heat-transfer  coefficients  are  larger  than  those  with  evapor¬ 
ation  at  a  body  surface,  all  other  things  being  equal.  A  decrease  in 
the  heat-transfer  coefficient  with  a  decrease  of  the  temperature  dif¬ 
ference  is  accounted  for  in  the  calculation  formulas  by  the  Gukhman 
number. 

It  was  shown  that  In  other  studies  [8,  9]  not  the  heat-transfer 
coefficient  a  was  calculated  but  the  over-all  heat-transfer  coeffi¬ 
cient  k  (see  formula  (39)),  which  decreases  with  an  increase  of  the 
temperature  difference.  Heat  and  mass  transfer  of  capillary-porous 
bodies  with  the  ambient  medium  is  a  single  interconnected  process  of 
heat  and  mass  transfer  In  the  layers  of  the  body  and  gas. 

DESIGNATIONS 

a  —  coefficient  of  thermal  diffuslvity  (m2/hr);  Cp  —  specific  iso- 
baric  heat  capacity  of  moist  air  (kcal/kg  •  deg);  ji  —  evaporation 
rate  for  flow  density  of  a  substance  (kg/ra2  •  hr);  1  —  a  character¬ 
istic  dimension  (m);  D  —  mass  (kg);  p  —  over-all  air  pressure  (mm  Hg); 
q  —  density  of  heat  flow  (kcal/m2  •  hr);  t  —  temperature  (°C);  T  — 
absolute  temperature  (°K);  r  —  specific  heat  of  evaporation  (kcal/kg);. 
w  —  velocity  of  air  motion  (m/sec);  a  —  heat-transfer  coefficient 
(kcal/m2  •  hr);  rj  -  coefficient  of  dynamic  viscosity  (kg/m  •  sec); 

X  —  coefficient  of  thermal  conductivity  (kcal/m  •  hr  •  deg);  v  —  co¬ 
efficient  of  kinematic  viscosity  (m2/sec);  p  —  density  (kg/m3); 
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Pio  —  relative  vapor  concentrations;  t  —  time;  i  —  distance  of  the 
evaporation  surface  from  the  body  surface  (m). 

Similarity  numbers:  Rex  —  local  Reynolds  number  value;  Pex  — 
local  Peclet  number  value;  Te-  —  local  Peclet  number  value  referred 
to  the  average  integral  velocity  wx  in  the  boundary  layer  fFex  = 

=  w  x/a);  Nu„  —  local  Nusselt  number;  Sc  —  Schmidt  number;  Pr  — 

Prandtl  number;  Gu  —  Gukhman  number. 

Indexes:  c  —  ambient  medium  (moist  air);  m  —  state  of  adiabatic 
saturation;  s  —  surface;  t  —  solid  body;  x  —  local  value  depending  on 
the  x-coordinate;  i  —  vapor;  2  —  dry  air;  0  —  heat  transfer  without 
mass  transfer. 
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